HEAT CONDUCTION AND ENTROPY PRODUCTION IN 
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Abstract. We investigate a class of anharmonic crystals in d dimen- 
sions, d > 1, coupled to both external and internal heat baths of the 
Ornstein-Uhlenbeck type. The external heat baths, applied at the bound- 
aries in the 1-direction, are at specified, unequal, temperatures Tl and 
Tr. The temperatures of the internal baths are determined in a self- 
consistent way by the requirement that there be no net energy exchange 
with the system in the non-equilibrium stationary state (NESS). We 
prove the existence of such a stationary self-consistent profile of tem- 
peratures for a finite system and show that it minimizes the entropy 
production to leading order in (Tl — Tr). In the NESS the heat conduc- 
tivity K is defined as the heat flux per unit area divided by the length 
of the system and (Tl — Tr). In the limit when the temperatures of 
the external reservoirs go to the same temperature T, k(T) is given by 
the Green-Kubo formula, evaluated in an equilibrium system coupled to 
reservoirs all having the temperature T. This At(T) remains bounded as 
the size of the system goes to infinity. We also show that the correspond- 
ing infinite system Green-Kubo formula yields a finite result. Stronger 
results are obtained under the assumption that the self-consistent profile 
remains bounded. 

Dedicated to Jiirg Frohlich and Tom Spencer with friendship and 

appreciation 

1. Introduction 

The rigorous derivation of Fourier's law of heat conduction for classi- 
cal systems with Hamiltonian bulk dynamics (or for quantum systems with 
Schrodinger evolution) with boundaries kept at different temperatures is an 
open problem in mathematical physics [9j. The situation is different for sys- 
tems with purely stochastic dynamics, e.g. for the Kipnis, Marchioro, Pre- 
sutti (KMP) model [H], where such results can be readily derived [T3l ITO]. 
An interesting area of current research are hybrid models in which the time 
evolution is governed by a combination of deterministic and stochastic dy- 
namics. The deterministic part of the dynamics is given by the usual Hamil- 
tonian evolution. The stochastic part can be of two different types. In the 
first type, the stochastic part is constructed to strictly conserve the energy, 
as studied in [5], or conserve also momentum, as in [l5 :2J. In the second type. 
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studied in [7J and [8j , the stochastic part is implemented by couphng the par- 
ticles of the system to "internal" heat baths with which they can exchange 
energy. To obtain a heat flow between external reservoirs at specified tem- 
peratures Tl,Tr, acting at the left and right boundaries of the system, the 
temperatures of the internal heat baths are chosen in a self-consistent man- 
ner by the requirement that in the nonequilibrium stationary state (NESS) 
there be no net energy flux between these baths and the system [71 Be- 
cause of this self-consistency condition, there is an average constant energy 
flux across the system in the NESS, supplied by the external reservoirs at 
specified, unequal, temperatures coupled to the boundaries of the system, 
and then carried by the Hamiltonian dynamics. A proof of Fourier's law 
for both types of hybrid models has been obtained for the case when the 
Hamiltonian dynamics is linear [5l[8], i.e., for a system of coupled harmonic 
oscillators. 

In the present work we investigate the self-consistent model for anhar- 
monic crystals. Unlike the case of the harmonic system, where it is known 
that Fourier's law does not hold when the "noise" is turned off (the heat 
conductivity then becoming infinite), one expects that in the anharmonic 
system with a pinning self-potential the conductivity will stay finite, i.e., it 
will satisfy Fourier's law, even when the strength of the noise goes to zero. 
We are quite far from proving this, however. What we do show here is that, 
for these anharmonic systems, conductivity for the finite system, defined by 
first letting both Tl and approach the same value, is given by a Green- 
Kubo formula. We also prove that this Green-Kubo conductivity is bounded 
in the system size, whenever the noise is finite. 

These results are obtained by studying the entropy production in the 
reservoirs in the NESS specified by the temperatures of all the reservoirs. 
We prove that the self-consistent profile minimizes, among all possible tem- 
perature profiles, the entropy production to the leading order in the differ- 
ence of the boundary temperatures Tl — T^. We then prove a uniform bound 
for the entropy production of a stationary state with a profile linear in the 
inverse temperatures. This leads to a bound on the leading term of the con- 
ductivity of the self-consistent system, given by the Green-Kubo formula for 
the finite system with all reservoirs at the same temperature T. 

Furthermore, we show that the corresponding Green-Kubo formula for 
the infinite system, giving the conductivity of the infinite system as a space- 
time integral of the energy-current correlations, is convergent. The bound we 
derive implies that the conductivity vanishes in the limit of infinitely strong 
coupling to the reservoirs. This behavior is also apparent in the explicit 
expression of the conductivity of the corresponding harmonic system (see 
equation (7.10) in [8j). The violent contact with the reservoirs most likely 
makes local equilibrium so strong that eventually no transmission is possible. 

There are no comparable results for anharmonic crystals with the first 
type of hybrid dynamics, but only some bounds on the conductivity [2|. 
Under the assumption that the self-consistent temperature profile remains 
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bounded, we show that the conductivity of the finite systems with a fixed 
Tl — Tr > is uniformly bounded in the size of the system. (This assumption 
is "clearly" correct but we are unable to prove it, see section 9.) 

The model considered is described in section [2] while section [3] contains a 
summary of the results proven in this paper. The existence of a NESS with a 
self-consistent temperature profile is proven in section HI Entropy production 
in the NESS is discussed in section \5\ and in section [6] we prove that the 
stationary state corresponding to the self-consistent profile minimizes, at the 
leading order in the temperature difference (Tl— Tr), the entropy production. 
Thermal conductivity in the NESS is discussed in section [7] and for the 
infinite homogeneous system in section [8l Finally, in Section [9] we present 
some concluding remarks. 



2. Time Evolution 

Atoms are labeled by x = {xi,...,Xd) G {-N,...,NY = An, N > 1. 
Each atom is in contact with a heat reservoir at temperature T^. The 
interactions with the reservoirs are modeled by Ornstein-Uhlenbeck pro- 
cesses at corresponding temperatures. The atoms have all the same mass 
m = 1. Their velocities are denoted by px and the "positions" by Qx, with 
Qx,Px S IS- We consider a mixture of fixed and periodic boundary condi- 
tions. The fixed boundary conditions are applied in the 1-direction, and the 
corresponding boundary sites will be used to make contact with external 
heat reservoirs. In the remaining directions, we apply periodic boundary 
conditions. Explicitly, let dAjy denote the set with |xi| = -|- 1 and let 
[x]- = —N + {xi + N) mod (2iV -f- 1), for i > 2. The boundary conditions 
are then qx = 0, for x G OAn- In addition, we let the inner boundary of Aat 
consist of those x with \xi\ = N, and we denote it by dA]y. 

As we will show, the heat flux in the stationary state will be entirely in 
the 1-direction and the properties of the system will be uniform in the d — 1 
periodic directions. We define A^ = {x G Aj^ : —N < xi < N} to label the 
bonds in the 1-direction. 

The Hamiltonian of the system is given by 



xgA 



N 



£x=§+t: ""^'^ ~ '^-'^^ \ ""^'^^'^ ~ + w{<ix), xEA 



N , 

(2.1) 



where the ej, i = 1, . . . denote the Cartesian basis vectors. We assume 
that V and W are smooth positive symmetric functions on M with quadratic 
growth at infinity: 

lim W"{±\) = VF^ > 0, lim V"{±\) = > . (2.2) 

A— >oo A— »oo 
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Clearly then, there are Ci, C2 > such that 

Ci{q^ - 1) < V{q) < C2{q^ + 1), Ci{q^ - 1) < W{q) < ^2(9' + 1) . (2.3) 

The dynamics is described by the following system of stochastic differen- 
tial equations: 

dqx =Pxdt, 

I (2.4) 

Apx = -dq^HN dt - jxPx dt + \l1-ixTx dwx{t) , 

with 7x > for all x £ A^r. Here 'Wx{t),x £ A^r, are independent standard 
Brownian motions (with average and diffusion equal to 1). The generator 
of this process has the form 

Ln= Yl i^p^'^Ndg, - dg^HNdpJ + Y 7x {Txdp^ - Pxdp,) 

x€Am xgAn (2-5) 

= ^ + 5, 

where A is the Hamiltonian part, anti-symmetric in the momentum vari- 
ables, and S is the symmetric part corresponding to the action of the reser- 
voirs. Then 

d 

Ln£x = ^ {jx-ei,x - jx,x+ei) + Jx, X £ An (2.6) 

i=l 

with Jx = 'JxiTx - pI) and 

Jx,x+ei 

= 0, if [x] Ajv or [x + Ci] Atv, (2.7) 
jx,x+e, = -\{P[x] + P[x+e,])y' {q[x+e,] " Otherwise. (2.8) 

In particular, then jx,x+ex can be non-zero only if x S A^. 

In section 3 of [16j it is shown that, for any choice of the temperatures T = 
{Tx > 0}, there exists an explicit Lyapunov function for the corresponding 
stochastic evolution, as long as 72; > for all x. This implies the existence 
of the corresponding stationary measure that we will denote by /u(T). 

If at least one Tx > 0, then the generator defined in ()2.5p is (weakly)- 
hypoelliptic, in the sense that the Lie algebra generated by the vector 
fields {A,dp^,x G Ajv} has full range in the tangent space of the phase 
space (M^'')'^^. In particular, the dynamics has probability transitions with 
smooth densities with respect to the Lebesgue measure on the phase space. 
If all Tx > 0, also the corresponding control problem has a strong solution 
(cf . section 3 in [16] , or [TT] ) and uniqueness of the stationary measure fol- 
lows from these properties. These methods could be extended to the case 
Tx > 0, at least if 'HN{p-,q) is strictly convex [l7j. The investigation of 
the uniqueness of the stationary measure goes beyond the purposes of the 
present paper, in particular, since zero temperatures will be relevant only 
in the general proof of existence of a self-consistent temperature profile in 
Section m So we will assume the uniqueness even in the case of temperatures 
not strictly positive. 
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The spatial periodicity will be exploited in the following to remove (most 
likely irrelevant) technical difficulties associated with irregular boundary 
behavior. To this end, we will assume that also the heat bath couplings 
respect this periodicity, i.e., we will always assume that ^y^ depends only 
on xi. Then in the case where also depends only on xi, the stochastic 
dynamics is fully invariant under periodic translations. Since the stationary 
measure /x(T) is unique, then also any of the corresponding expectation 
values must be invariant. 

We denote the constant temperature profile, Tx = To for all x G A^v, as 
Tq. Then /u(To) = ^xToi the Gibbs measure at temperature To, defined by 

/^To = exp(-WAr(p, q)/To)dpAq = Gt, {p, q)dpdq (2.9) 

We use as a reference measure and denote the related expectation by 

Oo- 

Computing the adjoint of with respect to the Lebesgue measure we 
hav^ll 

lT = -a+ y1 s:^'^ (2-10) 

xeAjv 

where Sx''^^ = 7^. [T^d^^ + l+p^dp^). We denote by /at = /iv(T) the den- 
sity of the stationary state /i(T) with respect to Lebesgue measure. This is 

the solution of /tv = 0. Due to hypoellipticity, is a smooth function 
of {p, q), and this implies also smoothness in T. To see this, note that dxy/w 
is the solution of the equation 

Lf^dTjN = -lydlfN. (2.11) 

Since the right hand side is smooth in (p, q) , this equation has a smooth 
solution, and smoothness in T follows by a standard iteration of the argu- 
ment. 

3. Summary of results 

Given the temperatures 0/j = {Qy}yi^R in a set i? C Atv, we say that 
a temperature profile T = {Tx}x(^An is self- consistent, if Tx = Ox for all 
x £ R, and the corresponding stationary state has the property 

lj)l)=Tx, for allxGAAr\i?, (3.1) 

where (•) denotes expectation with respect to the NESS, /u(T), assumed to 
be unique. Eventually we may choose R = dAj\f or part of it. But the 
following result is independent from the geometry. 

Theorem 1. For any choice of a non-empty R C A^r, and for any choice of 
temperatures 0/j = {Qy}y£R not all equal to 0, there exists a self- consistent 
temperature profile T = {Tx}xeAN\R- I""^ addition, if R and @r are invariant 

^We wish to reserve the standard notation for adjoint for certain weighted L'^-spaces, 
to be introduced later. Hence the notation L^^' for the adjoint here. 
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under translations in all of the d — 1 periodic directions of K^, then a self- 
consistent profile invariant under these translations can he found. 

The main body of our results concerns the case where the reservoirs on 
the two sides in the non-periodic direction are fixed to constant but unequal 
temperatures. We call this case the boundary layer setup. More explicitly, 
we then define R = dA^ = c^l-A-tv U c^rAtv, where SlAat = {x : xi = —N} 
and c^rAtv = {x : xi = A^}, and we fix on the left the temperatures Tx = Ti^ 
for X G SlAat, and on the right = Tr for x G SrAat, Tr < Tl. We also 
set /3l = -I and /Sr = Tp~^. Uniqueness of the self-consistent profile is not 
claimed in Theorem[Tl and this remains an open problem in the generality of 
the theorem. However, by restricting to small temperature differences and 
then relying on the implicit function theorem, we can get a self-consistent 
profile which is essentially unique. 

Theorem 2. For any given To > and N, there are eQ.,6 > with the 
following property: In the boundary layer setup with Ti^,Tfi such that |Tl — 
To\, |Tr — To| < there is a self- consistent extension of the temperature 
profile, T'^'^(Tl, Tr), and the extension is unique in the sense that no other 
profile T with maxa; \Tx — To\ < 6 is self- consistent. In addition, T*^^ is 
invariant under translations in all of the d — 1 periodic directions of A^, 
and the map (Tl,Tr) i— > T*'*^(Tl,Tr) is smooth. 

As an aside, let us remark that a careful inspection of the proof of Theorem 
[2] shows that its assumptions could be greatly relaxed, allowing for more 
general sets R and almost arbitrary potentials V and W. However, since 
the range of its applicability, determined by eq , can depend on and might 
go to zero as A^ oo, we have included the proof of the more general 
result in Theorem [TJ Furthermore, the assumptions about the asymptotic 
quadratic behavior of V and W will be used in latter proofs, and thus cannot 
be neglected. From now on, we assume that Tl — Tr is sufficiently small for 
applying Theorem [21 and let T'^'^ denote the corresponding self-consistent 
extension of the temperature profile, which is thus invariant under periodic 
translations and leads to a unique, periodically invariant, stationary state. 

For a generic profile T, we define the entropy production in a reservoir in 
the steady state //(T) as the energy fiux entering that reservoir divided by 
its temperature [IJ. The total steady state entropy production is then given 
by 

-(T)= ^ = E^^^^-i)- M 

By using the local energy conservation (j2.6p and denoting fix = T~^, we can 
write this as 

d 

^(T) = E E (f^-+-' - /5-)0x,x+e,) • (3.3) 

i=l xeA^ 

It is is well known [4j that a{T) > 0. 
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For the self-consistent profile T'^'^, there are no fluxes to the reservoirs for 
X OAn and consequently, as will be shown below, {jx,x+ei) = Jn for all 
X £ A'jy. The entropy production (13. 3p is then equal to 



a(T-) = (2iV + ir-'iP^ - P^)jN . (3.4) 

Thus we can estimate the magnitude of the self-consistent current by esti- 
mating the entropy production. 

Theorem 3. 

a{T'') < {2N + l)'^-2(/3a - /3l)'C^(T^', 7) (3-5) 
where, up to a constant c depending only on the potentials V and W , 

Consequently, the average self- consistent current is bounded by 

0<i^<C(T-,7) ^q^. (3.7) 

We expect, but are not able to prove, that the self-consistent profiles 
remain uniformly bounded in N. From such a bound it would follow that 
Jn = 0{N~^). We expect in fact that G [T^, Tl], as in the harmonic case 
[8j, c.f.. Section [9l What we can prove is that the first order term of Jn in an 
expansion in the imposed temperature gradient is 0{N~^). This is possible 
even without explicit knowledge about the asymptotics of the self-consistent 
profile. To this end, we consider also profiles T^^™ which are extensions in 
the boundary layer setup to a profile with linear f3x] we define 

(rf°)^' = ^^^^^i + /3R. + /?L), xgAn. (3.8) 
For these profiles, the entropy production satisfies 

= ^ {jx,x+e.)T^nn , (3.9) 



xgA'^ 



and we can derive a more precise bound for it. 

Theorem 4. Given b > 0, there exists a constant (^2(7; 6), depending only 
on 7, V, W , and b, such that for all < Tl < b, 

^(T/31-) < {2N + lf-\f3^-(3^fC2{r,b)- (3.10) 

Obviously, if To is any constant temperature profile, we have (t(To) = 0. 
Furthermore, ^-(Tq) = 0, and the second order derivatives can also be 
computed, yielding the following theorem. 
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Theorem 5. The Taylor expansion of a around a constant profile To at the 
second order gives 

c7(To + ev) = ^Q(v;ro) + 0(e3), Q(v;To)= Jy,ccVyV^ , (3.11) 



N 



where, with (•)o denoting the expectation in ;u(To), 

Jy,x = lx^y,x - lxly{hx{-LN{T^o)T^hy)Q, /l^ = -f - 1 . (3.12) 

The matrix J is positive, and if Jy^x is restricted to x,y E Ajv \ SAat, it 
becomes strictly positive. 

We now denote 5T = T^- Tr and To = (Tl + T^)/2. The next result says 
that the self-consistent profile minimizes entropy production, at least up to 
the leading order in the gradient of the imposed temperature difference, 5T. 

Theorem 6. The self- consistent profile T^^ is a smooth function ofTi^ and 
Tr. For a fixed To, its first order Taylor expansion 

T^^ = To + g'''ix)dT + Oi6T^), x£An, (3.13) 

is such that v = g'^'^ is the unique minimizer ofQ{'v; To) for fixed v{x) = 
X £ dAj\f, where we choose the +-sign for x £ SlAtv, and — for x G SrAat. 

Consequently, the self-consistent profile minimizes the entropy produc- 
tion up to errors of the order of 6T^. In particular, the leading term of the 
self-consistent profile can be obtained by minimization of the entropy pro- 
duction. This is consistent with the general belief that for small deviations 
from the equilibrium state imposed by external constraints, the stationary 
state will be such that it minimizes the entropy production with respect to 
variation in the unconstrained parameters [12] . The entropy production has 
also been studied by Bodineau and Lefevere [6] in this model, and originally 
by Maes, et al., [15] in the context of heat conduction networks. 

We define the thermal conductivity in the self-consistent stationary state 
(of the finite system) as 

2N+1 - 

K?5(ro)= hm — ^j^. (3.14) 



This is related to the entropy production by (|3.4p . yielding 

k'§{To) = Q(g^=;ro)/(2iV+i)rf-2 (315) 

where, as in Theorem [5l we have defined 

Q(v; To) = V J{To)w = T^ lim '"^'^"^^ . (3.16) 

Since g'^'^ minimizes Q(-), we find using (j3.10p 

^?^(To) < (2iV+l)2-'^r2 Jm <T-2C2(7;2ro). (3.17) 
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In particular, since the bound does not depend on N, this proves that the 
self-consistent conductivity defined by ()3.14p is uniformly bounded in N. It 
also has a Green-Kubo type of representation, as summarized in the follow- 
ing theorem. 

Theorem 7. The self- consistent conductivity is uniformly bounded in N 
and satisfies 

1 

K?^(ro) = -2 / Yl (-(2iV+l)Ve,g^^(x))0,,,.+ei(t)j0,e,(0))0 dt (3.18) 

° -eA^ 

where (•)o denotes the mean over the initial conditions distributed according 
to the equilibrium measure at the temperature Tq with the time evolution 
given by the dynamics corresponding to To, i.e., all the reservoirs are at 
temperature Tq. Here Veig^'^ix) = g^'^{x + ei) — g^'^{x) denotes a discrete 
gradient. 

A similar Green-Kubo formula can be obtained for the entropy production 
in the stationary state of the profile T'^^™. We will prove that 



xgA'j^ 



1 1 f°° 1 

1 + — )^y^ ^ (h,y+ei{t)j^,x+eM)o dt . (3.19) 



x,yeA' 



N 



By (I3.17p . this is always an upper bound for k^^{To). We expect the self- 
consistent profile to become linear away from the boundaries in the limit 
e — > 0, and to find Veig^'^(x) ~ — whenever xi is not too close to itA''. 
Although a proof of this property is still missing, we conjecture accordingly 
that both K^^{To) and the right hand side of (13.19P have the same limit as 
N ^ 00. 

The last result concerns the Green-Kubo representation of the conductiv- 
ity in the infinite system. Consider the infinite system on Z*^ with all 7a; = 7 
and all thermostats at temperature Tq. This infinite dynamics has a unique 
invariant measure given by the Gibbs measure on (M^'^)^ at temperature 
To, defined by the usual DLR relations. We denote also the infinite volume 
Gibbs measure by /xTo • The existence of the dynamics of this infinite system 
in equilibrium at any given temperature can be proven by standard tech- 
niques (cf. [18], where a similar result is proven for an analogous system in 
continuous space) . A proof of the existence of the dynamics in dimension 2 
for a certain set of non-equilibrium initial configurations is proven in pLOj. 
Consequently, we look at the dynamics starting from this equilibrium distri- 
bution, and let E denote the expectation over the corresponding stochastic 
process. 
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Theorem 8. There is a unique limit for 

-2 lim V / e-^*E [j.,x+ei(t)jo,ei(0)] dt = k(To) < - , (3.20) 

where C = E [(y(qei(0) — 90(0)))^] /To is finite and depends only on Tq. 

As we have mentioned in the introduction, the above bound for the con- 
ductivity goes to zero when 7 — > 00. 

As argued earlier, we expect the self-consistent conductivity and the 
Green-Kubo formula for the linear profile to have the same limit as ^ 00. 
However, inspecting the definition of the latter quantity in (j3.19p shows that 
this limit should be given by (j3.20p . provided the current-current correla- 
tions {jx,x+ei{t)jy,y+eiiO))o havc a Sufficiently fast uniform decay both in t 
and in the spatial separation |x — y| (the limiting infinite system dynamics 
are translation invariant also in the first direction, which should be em- 
ployed to cancel the sum over y in (j3.19p ). Therefore, we also conjecture 
that K^(To) k{To), at least along some subsequence of ^ 00. 

4. Self-consistent Profiles: Existence 

The following Lemma shows that zero temperatures cannot appear in self- 
consistent temperature profiles. (We will also give a second proof of local 
existence of self-consistent profiles in Section [6] which does not rely on the 
assumptions made about profiles containing zero temperatures.) 

Lemma 1. If{Tx,x G Aat} are not all identically zero, then (^Py'j > for 
all y G Atv- 

Proof. This is a consequence of the smoothness of the density of the transi- 
tion probability Pt{q' ,p'; q,p) of the process. Since / Pt{q',p'; q,p) dq dp = 1, 
for any {q',p') there exists an open set of positive Lebesgue measure A = 
A{Q\p'jt) such that 

/ Pt{q',p';q,p)dqdp>0. (4.1) 
J A 

If there exists x such that (p^) = 0, then 

= J fi{T-dq',dp') I plPtiq,p';q,p) dq dp (4.2) 
which clearly is in contradiction with (j4.ip . □ 

Proof of TheoremUi Given any collection of parameters u G [0, 00)^'^, x G 
R'^, let us define the corresponding temperature profile T(ii) by 

T{u)x = Tiu;&)x = h' '["^f' (4.3) 

I Ba;, it X e R. 

As before, we denote the density of the corresponding stationary measure by 
fN{q,P] T(ti), V, W). We have seen in the section [2] that, by the hypoelliptic 
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properties of the dynamics (cf. [16j), /n is a smooth function of {q,p) and 
consequently of T. By a straightforward scahng argument, we then have for 
any u and A > 0, 

X^'fNi^q, VAp; T(n), V, W) = fN{q,p; T(n)/A, Vx, Wx) (4.4) 

where Vx{q) = X^^V{'\/Xq) and Wx{q) = X^^W{^/Xq). An argument similar 
to that used at the end of section [5] to prove regularity in T shows that 
fNiq,P',T^iu)/X,Vx,Wx) is smooth in A. Under the conditions assumed on 
V and W, we have limA^oo ^a(q') = Vooiq) and limx^ooWx{q) = Woo(g) 
with VM = and Woo{q) = ^W^Lq^. 

We apply the scaling relation to prove that for high enough temperatures 
the system behaves essentially like a Gaussian. More precisely, consider 
arbitrary sequences A,i — > oo and b^") G [0,00)^^^, such that b^") converges 

to b G [0, 00)^^. Define further ri"^ = Xnb^x\ x G Kn- Then by the scaling 
relation (j4.4p . for any x', 

^(p2,)(TW,F,T^) = (p2,)(bW,y,„,ty,J {pl,){h,V^,W^). 

(4.5) 

The last expectation is with respect to the stationary state of a purely 
harmonic system. This system was studied in [8], where it was proved, in 
Sections 3 and 7, that there is a doubly stochastic matrix M, with strictly 
positive entries, such that for any profile of temperatures b and for all x' , 

{pI.) {h,Voo,Woo) = Yl M^'yby 

(Strictly speaking, the result was proven only for periodic profiles in [8]. 
However, the above properties, linearity in b, as well as positivity and double 
stochasticity of M, are easily generalized for non-periodic profiles, although 
we do not go into details here.) Since M^y = 1 for all x, this implies 

(pI,) {h,Voo,Woc) < maxby = \\h\\^ , (4.6) 

and the equality holds if and only if b is a constant vector, i.e., bx is inde- 
pendent of x. 

We can now prove the existence of a self-consistent profile. Let i?^ = 
Ajv \ R, and consider the mapping F : X ^ X, X = [0, 00)^' defined for 
u G [0,00)'^'', X G R^, by 

Fiu), = {pl)iT{u),V,W). (4.7) 

Since some of the temperatures are kept fixed to non-zero values, the hy- 
poelliptic properties of imply that F is everywhere continuous. For 

any L > define Xl = [0, L]-^" C X. We will soon prove that there is an 
L > such that F{Xl) C Xl. Since Xl is homeomorphic to the unit ball 
of M}^"^ and F is continuous on Xi, we can conclude from the Brouwer fixed 
point theorem that there is at least one u G X^ such that F(u) = u. By 
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LemmalU if there is x such that = 0, then F{u)x > 0, and such u cannot 
be fixed points. Thus for any fixed point < < < c)o for all x, and 
T(u) is then a proper self-consistent temperature profile. 

We prove the existence of a constant L, for which F{Xl) C Xi, by 
contradiction. If no such L exists, then for all L > there is u^^-* G Xl 
such that ||F(n(^))||oo > L- Then necessarily ||ti^^^||oo — > cO; since otherwise 
there would exists a convergent subsequence, which is incompatible with 
\\F{u'^^^)\\oo oo. Let \l = llii^^^lloo and = \l^u^^\ so that Al ^ oo 
and 1 1 u^'^-' I loo = 1- The sequence {v^^^) belongs to a compact subset of X, 
and we can find a subsequence such that v^^^ v m. X. For this final 
subsequence we can apply (|4.5p and (j4.6p . which shows that for all x 

limsupA;^^F(ALi;(^))^ < ||i;||oo = 1. (4.8) 

L 

Equality is not possible here, as the limit b of Xj^^T{\iv^^'') has at least 
one component equal to one, but bx = for all x £ R, and thus b cannot 
be a constant vector. However, by construction, for every L there is x{L) 
such that F{\lv'^^^)x[l) > L > ||ti(L)||oo = A^, which leads to contradiction. 
This proves the existence of L > with the required properties and concludes 
the proof of the first part of the theorem. 

For the second part, let us first point out that, if R is invariant under 
all periodic translations of A^v, it must be of the form R = Ri x 
where In = {—N,---,^} and Ri C In is non-empty. Similarly, @x can 
only depend on xi. Let R^ = In \ Ri, let Pi denote the projection on 
the first axis in Z^, and define R' = PiR'^ = R\ x {0}, which is a subset 
oi R'^ = Kn \ R. If -R' is empty, R = Kn and there is nothing to prove. 
Otherwise, let us consider the map F' -.X' ^ X', X' = [0, oo)-^', defined by 
F'(n), = (p2)(T'(n),F,Vr), where 

r».. = ^^-' (4.9) 
I Qx, otherwise. 

Every such T'{u) is clearly invariant under all periodic translations. We can 
then repeat the analysis made above for F' and conclude that it has a fixed 
point u with < < oo. Since T = T'{u) is periodic, the dynamics is 
completely invariant under periodic translations, implying that also expec- 
tation values in the unique stationary state are invariant. Therefore, for any 
X G R'^, we have {px){T) = (pp^^)(r) = up-^x = Tx- This proves that T is an 
invariant, self-consistent profile. □ 



5. Entropy Production Bound 

In this section we prove the entropy production bounds stated in The- 
orems [3] and [H Given a generic profile of temperatures T, we recall the 
notation Jn = /7v(T) for the density of the stationary measure /i(T) with 
respect to Lebesgue measure, and let (•) denote expectation with respect 
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to /u(T). A simple computation shows that {Alnf]\f) = ioi A defined in 
()2.5p . Therefore, by stationarity we have 

= - (Liv In /at) = - ^ (5a: In /at) (5.1) 

X 

where = '^x{Txdl^ -PxdpJ. Let = In/Gt^, where Gt = Z^^e"'^'^/^, 
as in (12.91). Then we can rewrite the last term as 



-{SxlnfN) = -j {Sxlnil!x)tpxGT^dpdq- j SxilnGrJfNdpdq. (5.2) 

Since PxGt^ = -Txdp^Gx^ and SxihiGrJ = -^x{Tx - pl)/Tx = -Jx/Tx, 
we find by integration by parts that 

-{SxlnM=Txjx [^^^^GTApdq+^. (5-3) 
J Wx Tx 

So by (j5.ip . the entropy production satisfies 



2?. = 7xr. / i%M!G^^dpdg. (5.5) 



where 

Ipx 

In particular, (t(T) > 0, and by using the local conservation of energy, (j2.6p . 
(I33D holds. 

Let us for the remainder of this section assume that T is a temperature 
profile which is invariant under the periodic translations. The results then 
hold for both and T'^"'^. Obviously, then by 1^ 

a(T) = ^ {Px+e, - Px){jx,x+e,) (5.6) 
xGA^ 

Therefore, it will suffice to find a bound for \{jx,x+ei)\- 

Applying the definition of the current observable, (|2.7p and (|2.8p . and 
then integration by parts, shows that 

1 f ^ 

{jx,x+ei) = -^ / V'irx)YiJx'Px'GT^,\^,^^^^^^dpdq 

n=0 

= -E^ / y'{rx)GT^,dp^,i^x'dpdq (5.7) 

Z J x'=x+nei 

n=0 

where rx = qx+ei — Qx- We use that 1 = 'tpx^'^ /il^xC^ whenever i/^x' / 0, and 
then apply the Schwarz inequality. This shows that 

\{jx,x+e,)\^ < max ^ {V'irxf) \Y,Vx+ne, . (5.8) 



n=0 
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Therefore, we have obtained the following relation between the total sum of 
currents and the entropy production 

(E IO-...-fe.)l)'<max| E {V'ir.f) ^ ^^^^-^^ 
xgA'^ " xeA'j^ xeA'j^ n=0 

< a{T) max ^ V {V'{r,f) . (5.9) 

^^^^ "^y .t^. 

For this bound to be useful, we still need to consider J2xeA'j^ {'''x)'^) ■ 
Since LN^q'^) = 2qxPx, we have {qxPx) = for all x. Similarly, = 
^xeAN^xiTx-pl) implies Ex Tx^x = Ex7x'(Px)- Now 

LNi E ^'^'^^O = J2 P^^~ Yl l^^i^'^- E TxPxgx, (5.10) 



and thus 



E ^^^^ ^ mm7y E (^^x) = min7j^ / ^ <lxdq.'H) ■ (5.11) 

xgAiv xgAiv \xeAjv / 

From the asymptotics of V and we can conclude that there are C > 
and C" > such that 

V'{rf < C{rV'{r) + C) and rTy'(r) > -C" . (5.12) 

But since 

d 

d,^n = W'iQx) + E {y'id- - l-e,) - y'iQ-+e, - qx)) 

i=i 

+ ^ (1(:e G d^AN)V'{-qx) - l(x G aLAjv)l/'(gx)) , (5.13) 
with 1 denoting the characteristic function, we have 

d 

J2 Qxd,^n= J2 [lxW'{qx)+Y,^^'i^ 

xeA'^ xS^rA xg9lA 

> ^ r,y'(rx)-|Aiv|C'((i + l). (5.14) 

xgA'j^ 

Combining this with (j5.1ip shows that 

E {y'irx)') < c\AM (c'{d + 2) + "^^""y^y^y ) . (5.15) 
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Consequently, there is c > 0, which depends only on V and W, such that 

Let us next consider the case T = T^'™. Applying the definition of T^'™ 
to dSS]) shows that then 1^ holds, i.e., a(T^i™) = ExeA^ 0'^',a:+e,)- 

Then by (I5.16P and using the fact that T^^™ < Tl 

aiT"''^) < c' - (2iV + + TO' , (5.17) 

where c' is a constant depending only on 7, 1/, and M/^. Therefore, we have 
now proven Theorem [H 

Finally, let us consider the self-consistent case, T = T'^'^. For the corre- 
sponding stationary measure we find from (|2.6p , 

d 

{{j^,x+e,) - {jx-e,,x)) =0, X Mat . (5.18) 

i=i 

Since the system, including the self-consistent profile, is periodic in any of 
the Cartesian directions Cj, i > 2, also the unique stationary measures are 
invariant under translations in these directions. Therefore, 

{jx,x+e,) = {jx-e„x) , i / 1, X G Aat . (5.19) 

Consequently, by (I5.18P and (12. 6p . 

{jx,x+ei) = {jx-ei,x) , X ^ OAn , 

{jx,x+ei) = (Jx) = ixin - (pI)), X e d^AN , (5.20) 

{jx-ei,x) = -{Jx) = 1x{{pI) - T^), X £ SrAat . 

We denote the constant current by Jn, i-e., now we have {jx,x+ei) = JN-, for 
all X S A^. Therefore, by (15. 6|) . 

a(T^^) =3nY1 - (^-) = - /^0(2A^ + 1)'-' , (5.21) 

xgA'j^ 

which proves (j3.4p . This immediately implies that sign(TL — Tj^)jN > 0. 
But on the other hand, Jn = jx^YlxGA'^U^,^+ei)-' ^^^^^ ^^^^ ^'^^ 
self-consistent profile a{T^^) = ExeA'J^^,^+e^) ■ Applying dESI) then 

completes the proof of Theorem [3j 

6. Minimization of entropy production 

For a given To > 0, we use the Gibbs measure = GT^dpdq as a refer- 
ence measure and we denote the related expectation by (•)q. We consider the 
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generator L on the Hilbert space L?{iiTa)- Recall that for any temperature 
profile T = {T^^x G Aat} we have L = L(T) = A + S{T). Its adjoint is 

L* = -A+Y,Sl (6.1) 

xeAjv 

where 5^; = 7^; [T^d^^ - Pxdp^), and thus 

AT,., 

with ATx = Tx — To and 



S; = Sx + i-^Pxdp^ + hx) (6.2) 



hx = ^-l. (6.3) 

Observe that {hxhxi)o = '^^x,x' and —Sx,Tohx = '^Ixhx- 

Set Lo = L(To) and consequently = + S'(To)* = -A + S'(To). 

Lemma 2. For all y, x 

'9t^(p^)mt)It=To = ly{hyi-^o)~^^^)o = 7y{hx{-Lo)~^hy)o . (6.4) 

Proof. Let us denote by / = /(T) the density of fJ,{T) with respect to ;UTo. 
Then / is solution of the equation L*(T)/(T) = 0. Since the coefficients 
in L*(T) are smooth in T, / is smooth in T and fy = dTyf{T^) solves the 
equation 

L*(T)/,(T) = -idT,L*)iT)fiT) = (Todl - 2p,9,, + hy) /(T) . 

(6.5) 



Since /(To) = 1, we have found that /^(To) is solution of 

-L*fy{To) = ^hy. (6.6) 

Notice that fy has a bounded L^^fiTa) norm (cf. [20]), and by a standard 
argument (multiply equation (j6.6p by fy and integrate with respect to /i-To) 
we obtain a bound 

Y.^x{{dpJyf)o<iyT-\ (6.7) 

X 

Now, since KGt^ = -dp^p^Gr^), 

{hx)^iT) = {Pxdp J {T))o. (6.8) 
Then differentiating with respect to Ty we have 

dTy{hx)^(T) = {PxdpJy{T))o (6.9) 
and taking the limit T Tq 

9t,(/ix)mt)|t=t. = MJy{To))o = {hxfy{To))o = ^Jhx{- L*)'' hy) o 

(6.10) 

Observe that, since h is an even function of p, one can, by a change of 
variables p — > —p, replace Lq with Lq in (|6.10p . This proves (16. 4p . □ 
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Define F : Wl"" M^^ as 

F,{T) = (J,)^(T) = 7. {Tx - {pI)^(t)) ■ (6.11) 
Its Jacobian at T = To is given by 

Jy,x = lx5y,x-lx (Px)mT) I T=To = Ix^y^x- lxly{hx{- LqY^ hy) q. (6.12) 

Observe that J is symmetric and that -F(To) = for any value of Tq. It 
follows that is an eigenvalue of J , and we will show shortly that i7 > 0, 
and the eigenspace corresponding to is one-dimensional and generated by 
the constant vector. Then the matrix M = {Jx,y)x,yeR'= is invertible, and 
thus there is a neighborhood in (Tl,Tr) containing (To, To) such that the 
implicit function theorem can be applied to obtain a self-consistent profile. 
This implies that constants Eq and 5 for the first part of Theorem [2] can 
be found. It also follows that T'^'^(Tl, Tr) is smooth. To see that it must 
also be invariant under the periodic translations, we first point out that in 
the boundary layer setup clearly any translate of a self-consistent profile 
is also self-consistent. Since the translations correspond to a permutation 
of indices, they remain in the neighborhood determined by 5, and thus by 
the uniqueness of the self-consistent profile in this neighborhood, T'^'^(Tl, T^) 
must itself be invariant. 

Therefore, to complete the proof of Theorem [2] we only need to prove the 
following Lemma. 

Lemma 3. ^7 > 0, and J a = implies is a constant in x. 

Proof. Let a E M^'^, and define h = '^^^Km follows from the 

antisymmetry of A and the symmetry of Sq: 

{iAh)i-Lo)-\Ah))o = {hi-So)h)o - {{Soh)i-Lo)-HSoh))o . (6.13) 

Since Soh = —2Y^^ax^xhx, we obtain 

{{Ah){-Lo)-\Ah))o = (/i(-5o)/i)o -4^a^ay7^7y(/i^(-Lo)"^/iy)o 

x,y 

= 4^7j,.a^. - 4:'^axay-/x-fy{hx{-Lo)~^hy)o 

X x,y 

= i'^axttyJx^y . (6-14) 

x,y 

Therefore, to prove that J has the properties stated above, it suffices to 
study the left hand side of (j6.13p . and to prove that it is always positive, 
and equal to zero if and only if a is a constant vector. (Studying real vectors 
a suffices here, as ^7 is a symmetric matrix.) 

In fact, define u = {—LQ)~^{Ah). Since for any observable F belonging 
to the domain of A, {F{AF))q = 0, we have then 

{iAh)i-Lo)-\Ah))o = {u{-So)u)o = Y.lxTo{{dp^uf)o > 0. (6.15) 
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This proves the required positivity. In addition, if the left hand side is zero, 
then u{p, q) cannot depend on p, and thus 

- Lqu = -Au = - '^Pxdg^u{q) = Ah = - —'^ o.xPxdg^7i (6.16) 



It fohows, for ah x. 



Thus the function 



-a^dg^n = dgMQ) ■ (6.17) 



2 

To 



QiQ) = Y«(9) - E ^-W(l-) (6-18) 

xgAn 

satisfies, by (|5.13p . 

d 

dq^Siq) = [E (^'(^a; - Qx-e,) - V'iq^+e^ - q^)) 
i=i 

+ ^l(x e d^AN)V'{-q,) - ^l{x e d^AN)V'{q,)\ . (6.19) 

For X G and A; = 1, 2, ... we differentiate (I6.19j) with respect to qx+e^ 
and obtain 

- axV"{qx+et: - qx) = 5g^,g,+,^^(9) = -ax+ekV" {Qx+eu - Qx) ■ (6.20) 
Since there exists an rg such that V"{rQ) > 0, this imphes a = const. □ 

We can now conclude that for any Tq > 0, there is eo > such that for all 
|e| < eo a self-consistent profile corresponding to Tl = To+|, Tr = Tq — | can 
be found. This profile is differentiable with respect to e and the derivative 
satisfies for x dA^ 

= ^/x{T{e; To)) = J] ^5T,F.(T(e; T,)) . (6.21) 
Therefore, we have J2yeAjs, •^x,y^^-^^ = 0. This shows that for x dA^, 



dT,{e;To 



de 



= {M-%/-( ^y,y'- E ^vA^ (6. 

y09Ajv \'GORA^r y'ddi^hN ^ 



where M = {Jx,y)x,y^dAjs, is a strictly positive matrix, and thus invertible. 

Recall the definition of entropy production given in (13. 2p . By (j5.4p we 
have then always cr(T) > 0, with equality when T = To, a constant profile 
given by To > 0. Since 

'^''^x J^x 
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we have for the constant profile 

da , d 



^^(T„) =-T-S. + T-.^(X:7„(p?>)^^^_. (0.24) 

As mentioned earlier, for any profile Ylyly {Py) — Yly'ly'^y^ thus we 
have proven that 

^(To) = 0. (6.25) 
A similar, but a slightly longer calculation, shows that 

dT dT ^'^"^ ~ j^^'-^^'y ~^ "^V'^) ~ '^'^^^y ■ (6.26) 

By dividing Atv into i? 7^ (the fixed thermostats) and R'^, we can con- 
clude from the previous results that the symmetric matrix M = {Jx,y)x,y&R'' 
is strictly positive. By ()6.25p and (I6.26p . the Taylor expansion of a around 
To yields 

a(To + ev) = ^ Yl ^x,yVxVy + 0{e^) . (6.27) 

° x,y£AM 

This proves Theorem [5l For fixed e and Vx, x £ R, the quadratic form 
corresponding to the leading term has a unique minimizer, given by 'ui™™'* = 
Vx, X £ R, and 

^(min) = - ^ ^ {M~\yJy,y>Vy> , ioi X £ R'' . (6.28) 

y&R" y'eR 

Let us next consider the case studied earlier, with the opposite boundaries 
fixed at two different temperatures Tl and Tr. Denote 6T = Tl — Tr, which 
we assume to be positive, and To = (Tl + T^)/2. Let us consider a sequence 
of Il,Tr for which Tq remains fixed and 6T — > 0. We assume that T is 
a sequence of profiles with boundary values on R equal to Tl and Tr, and 
which has a Taylor expansion 

Tx = To + g{x)5T + 0{5T^) (6.29) 

where g'lsa, function for which g{x) = 1/2 for x £ di^Aj^ and g{x) = —1/2 for 
X G c^RAiv. By (j6.27p . the entropy production will be of the order ((5T)^, and 
the leading term is minimized by g{x) = v^^^^ corresponding to Vx = 
with +, if X G Oi^An, and — , if x E 5rAjv. 

We have proven in the beginning of this section, that the self-consistent 
profile can be chosen for all sufficiently small 5T so that it is differentiable in 
the boundary temperatures. In particular, comparing (16.221) to (j6.28p shows 
that 

T^^ = To + g'^{x)5T + 0{6T^) (6.30) 
with g'^'^ = v^™™\ We have thus proven Theorem [6l 
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7. Conductivity of the Finite System 

In the following we again set A'^ = AjvVSrAat, and consider, as in Section 
[5l a generic profile T which is invariant under periodic translations. Let (•) 
be the expectation with respect to the corresponding stationary state. It 
is convenient now to use as a reference measure the inhomogeneous Gibbs 
measure vt = G{T; q,p)dqdp, with 

GiT-q,p) = '-^^^tZ^4^MlIlA (7.1) 

Zj 



where Ex is defined in (j2.ip . Notice that S is automatically symmetric with 
respect to i't, while the adjoint of A is given by 

[l^ ;^)jx,x+ei . (7.2) 

JV 



Let us next inspect T = T^"^ and denote by / the density of the self- 
consistent stationary state with respect to v^^si. Let us fix To = "^^^^ 
with e = = Tl — Tr > 0, as before. Repeating the argument used in 
section [2l we find that / is smooth in e, so a first order development in e 
is justified. Using the expansion of the self-consistent profile, (|6.29p . shows 
that u = d^f\^=Q is solution of the equation 

{-A + 5(To))n = ^^C^i.,.+ei . (7.3) 



Explicit formulae for the derivatives of the self-consistent profile, g'^'^(x), are 
given in ()6.22p . 

Recall the definition of the conductivity of the finite system, ()3.14p . Since 
we have already proven Theorems l-[6l the argument given before Theorem 
[7] in Section [3] provides a proof that Kj\f{To) is bounded in A^. On the other 
hand, by (f7^ . 

^7v(To) = hm ^^"^ ^ ^ I Yl 



ATI 



xGA^ 



= J™n ^^0.^^) = + l)(^.?0,ei>0 • (7.4) 

Define u{q,p) = u{q,—p), and observe that, since jx,x+e-i is antisymmetric 
inp, 

{A + S(To))n = -Y ^^C^J---+e, . (7.5) 



xGA'^ 
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Thus 

POO 

HNiTo) = -i2N + l){ujo,e,)o = {2N + l) / dt{u{t)jo,e,mo dt 

Jo 

1 

= r^l^ (-(2^ + l)Ve,g^^(x))(j.,.+e,(t)jo,e,(0))o dt 



(7.6) 



xgA^ 



where (•)o denotes taking the initial data distribution according to the 
equihbrium measure at the specified temperature To, and then considering 
the time-evolution corresponding to the stochastic process with all heat- 
bath temperatures set to Tq. We have used here the property that then 
(^(i)jo,ei (0))o — > ('u)o(jo,ei)o = for t ^ cxD. This completes the proof of 
Theorem [71 

Repeating the same steps for T = T^'™, for which deTx^^^\^^Q ~ ~§n^ 
proves also the validity of (j3.19p . 

8. Conductivity of the Infinite System 

We prove here Theorem [8] concerning the infinite system on (M?'^)^'' with 
all 7x = 7 and all thermostats at temperature Tq. This infinite dynamics 
has a unique invariant measure given by the Gibbs measure on (M?'^)'^'' at 
temperature To, defined by the usual DLR relations. We denote this measure 
by iJ,Ta and its expectation by (•)q. Consequently we look at the dynamics 
starting from this equilibrium distribution. 

We adapt here an argument used in [3]. Introduce on LP'{^To) a degenerate 
scalar product 

i^))=Y. [{w.i^)o - {v)o Wo] > (8-1) 

where Tx is the translation operator. The scalar product can also be obtained 
via the limit 

= lim Cov^ ($„(^,$„V) = lim {{^nf'^ni^)o - (^>„¥j)o(^nV'>o) , 

n—>oo " n—>oo 

(8.2) 

where maps functions into the corresponding "fluctuation averages" in 
A„,, a square box of linear size n centered at 0. Explicitly, 

i'^ni^){q,p) = (^-V')(?,P) • (8.3) 

V x€A„ 

The scalar product is degenerate, since every function of the form (f> = 
ip — Txip is in its kernel. We denote by the corresponding Hilbert space 
of square integrable functions. More precisely, is a space of classes of 
functions such that each of its elements can be identified with a function in 
L^(^To) up to a translation. 

Observe that A and S are still respectively anti-symmetric and symmetric 
with respect to the scalar product ((•,•))• We also introduce the semi-norm 

Ml = {{^,{-S)y,)) (8.4) 
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and let Hi denote the corresponding Hilbert space obtained by closing 
with respect to || • ||i. To see that is a semi-norm, in particular, that 

it is positive, we can employ the easily derived identity 

\\ip\\l= lim {^r.vi-S^n^))o . (8.5) 

Since 5 acts only on velocities, || • ||i has a kernel consisting of all functions 
which depend only on q, the position variables. Thus also TCi is a space of 
equivalence classes of functions. 

Let A > be given and let ux be the solution of the resolvent equation 

Xux - Lux = JO.ei • (8.6) 

The solution can be given explicitly in terms of the the semigroup P* gen- 
erated by L = ^ + 5, 

/•oo 

uxiq,p)= e-^\P'jo,e,)iq,p)dt. (8.7) 
Jo 

Obviously, 

Co := ((jO.ei, JO,ei)) = X] {jo,eJx,x+ei)o < To {{V {qe^ - Qo))^)o < ^ > 

(8.8) 

and thus jo,ei £ Then ux G L^(/xj'^), and by stationarity (ma)o = 0. We 
will show next that, in fact, ux G TCi- From (18. 6p we obtain 

X{i^nUxf)o + {{<^>nUx)i-S){'^nUx))o = {i'^nUx)i<^njo,e,))o , (8.9) 

where we have used translation invariance of L and antisymmetry of A. 
Since Sjo^ei = — 7jo,ei , an application of Schwarz inequality yields 

(($„UA)(^„j0,ei))0 = 7"^((^nMA)(-5)(^ni0,ei))0 

<^-\{^^Ux)i-S)i<^nUx))l^\i^nj0,e,){-S)i^nj0,eX^^ 
= T'^'{i<^nJO,e,)X^\{<^nUx)i-S)i<^>nUx))l^' . 

(8.10) 

Consequently, we have 

{i'^nUx){-S){<^>nUx))o < 7"' ((^nj0,e j')o ™ l~'Co , (8.11) 

which implies that 

X{{ux,ux)) <-f-^Co (8.12) 

and 

llnAll? <7"'C'o. (8.13) 

Therefore, ux G Tii and by (I8.13p . we can extract a subsequence, which we 
still denote with ux, weakly convergent in TCi to uq. 

Let ux{p,q) = u\{p,q) + u'^{p,q) where u\ and are respectively sym- 
metric and antisymmetric in the p's. Since jo,ei is antisymmetric in the p's, 
we have that {{ux, jo^ei)) = {{''^'x^ h,ei)) ■ Furthermore, 5 preserves the parity 
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in while it is inverted by A. So we can decompose the resolvent equation 
as 

Xul - Sui -Au'{ = 0, 

- 5< - Aul = jo,ei . ^^'^^^ 

Taking a scalar product of the first equation with u^, of the second with u^, 
and using the antisymmetry of A, we find 

(«, jo,e,)) = HK,ul)) + {{Ul, i-S)K)) - {{u'l,Aut)) 

= ^(«,<)) + A(K,<)) + ((u,,(-5K)). ^ 

Since 

^ ul{p,q)fiTAdp)=0, (8.16) 

where ft-Toidp) is the centered Gaussian product measure of variance To, and 
S has a spectral gap 7 in L'^{flTo{dp)), we have that 

(K,^)) < -{{u^,{-S)ux)) < Cor'- (8.17) 

7 

In particular, Uq £ C? . Thus by taking first the limit as A —> we have 
X{{u^^,ul)) 0, then as ^ we have v{{u'^^Uq)) 0, and finally we 
obtain from ()8.15p 

((^tO,JO,ei)) = {{uq,{-S)uo)) = \\uq\\1. (8.18) 
On the other hand, we have 
((^io,jo,ei)) = lim((nA,jo,ei)) = hm [A((nA,nA)) + ((ua, (-S')nA))] 

A- — ^0 A — 



> lim A((nA,UA)) + \\uo\\\ 

A- — >0 



^.19) 



This implies 
and 



limA((nA,nA)) = 0, (8.20) 

A- — *0 



ll^iAlli Ikolli ■ (8.21) 
Therefore, u\ — > uq strongly in "Hi. 

Uniqueness of the limit follows by the following standard argument. Sup- 
pose that A„ is the chosen subsequence such that converges to uq, and 
suppose Vm is another sequence such that u^^ converges to uq. Then, simi- 
larly as we have done in equation ()8.15p 

(K„,io,ei)) = z^m(«„,<J) + A„((7xi^,<J) + ((UA„,(-S)n,^)) (8.22) 
which implies 

(«,jo,ei)) = ((txo,(-5)uo». (8.23) 
Using n"^ instead of ti^^ , we find similarly 

{{ul,jo,e,)) = {{uo.{-S)uo)). (8.24) 

Combining these with (|8.19p shows that ||uo — ttolP = 0, i.e., uq = uq. 
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Thus the conductivity k,(To) defined by (I3.20p is independent of the sub- 
sequence chosen for A. Moreover, we have 

<T^o)=T^ ((^io,Jo,ei)) = ll-uolli < 7^ < 7^ (8.25) 

This completes the proof of Theorem [HI 



9. Concluding remarks 

While all the results obtained in this paper are as expected, the difficulty 
of actually proving things about the NESS of systems with nonlinear dy- 
namics is immense. This is well illustrated by the impossibility (for us) of 
obtaining a bound on the self-consistent temperature T of the second oscil- 
lator in a system consisting of three oscillators with T\ = Tl, T3 = T^, and 
the Hamiltonian is as in (|2.ip with 7a; = 7 > 0. We certainly expect that T 
will satisfy < T < Tl, but do not know how to prove this. All we know 
is that there exists a T = and that j = Tl — (p^) = — Tr > 0. 

We also know for general that when Tl,Tr To, then there is a self- 
consistent choice T To, and that this in this limit (2iV + 1)j/(Tl — Tr) 
is bounded and given by the Green-Kubo formula ()3.18p . Beyond this how- 
ever we are stymied except when V and W are harmonic. In that case T is 
given by (I3.13P without a correction term for any Tl, Tr, and due to explicit 
expressions g'"^(x) can be analyzed in great detail, proving Tr < T < Tl. 
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